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Abstract. Chronic wounds represent a major public health problem affecting 6.5 million people in the United 
States. Ischemia represents a serious complicating factor in wound healing. In this paper we analyze a recently 
developed mathematical model of ischemic dermal wounds. The model consists of a coupled system of partial 
differential equations in the partially healed region, with the wound boundary as a free boundary. The extracellular 
matrix (ECM) is assumed to be viscoelastic, and the free boundary moves with the velocity of the ECM at the 
boundary of the open wound. The model equations involve the concentrations of oxygen, cytokines, and the densities 
of several types of cells. The ischemic level is represented by a parameter which appears in the boundary conditions, 
< 7 < 1; 7 near 1 corresponds to extreme ischemia and 7 = corresponds to normal non-ischemic conditions. 
We establish global existence and uniqueness of the free boundary problem and study the dependence of the free 
boundary on 7. 
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1. Introduction. Wound healing represents the outcome of a large number of interre- 
lated biological events that are orchestrated over a temporal sequence in response to injury 
and its microenvironment. The process involves interactions among different soluble chem- 
ical mediators, different types of cells, and the extracellular matrix (ECM). Among the var- 
ious factors that affect the healing of a wound, the tissue oxygen level is a key determinant 
lfTTll26l . Although hypoxia is generally recognized as a physiological cue to induce angio- 
genesis |U [25] EH [14), severe hypoxia cannot sustain the growth of functional blood vessels 

G2HUIIS1III1I221. 

There have been several mathematical models of wound healing which incorporated the 
effect of angiogenesis ll20l [191 l3l l24l . Mathematical models of angiogenic networks, such as 
through the induction of vascular networks by vascular endothelial growth factors (VEGFs) 
J5] [6), were developed by McDougall and coworkers |[T6l [27], based in part on the work 
of Anderson and Chaplain J2|, in connection with chemotherapeutic strategies. The role of 
oxygen in wound healing was explicitly incorporated in the works of Byrne et al. (3 1 an d 
Schugart et al. Il24l . In particular, it was demonstrated in [24] that enhanced healing can 
be achieved by moderate hyperoxic treatments. In ll22l . the impairment of dermal wound 
healing due to ischemic conditions was addressed in a pre-clinical experimental model. In a 
more recent work 11281 . Xue, Friedman and Sen developed a mathematical model of ischemic 
dermal wound-healing. The model consists of a system of PDEs in the partially healed region 
which is modeled as a viscoelastic medium with a free boundary surrounding the open wound. 
Simulations of the model were shown to be in agreement with the experimental results in [22 1. 

In this paper we study the model in 11281 by mathematical analysis. In particular we 
prove that the free boundary problem developed in that model has a unique global solution, 
and that the open wound does not close under extreme ischemic conditions. We also show, by 
simulations, that non-ischemic wounds do heal. In Section[2]we formulate the mathematical 
model for a radially symmetric geometry as in ll28l . The ischemic level is determined by a 
parameter 7, < 7 < 1; 7 near 1 corresponds to extreme ischemia and 7 = corresponds 
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to normal non-ischemic conditions. In Section[3]we show that the free boundary is monotone 
decreasing, and in Section|4]we derive a priori estimates. In Section[5]we transform the free 
boundary problem into a problem in a fixed domain; this is a convenient form for proving, 
in Section [6] local existence and uniqueness of a solution. The extension of the solution to 
all t > is also established in Section [6] by using the a priori estimates derived in Section 
|U In Section [7] we consider the case of extreme ischemia (namely, 7 near 1) and prove that 
the wound's boundary stops decreasing after some finite time. In Section [S] we establish 
some properties of the solution for wounds that do not heal. Section [9] simulates the radius 
of the wound when the parameter of the system are chosen, as in ll28l . based on biological 
literature. The simulations suggest the following conjecture: there exists a parameter 7* such 
that wounds heal if < 7 < 7* and do not heal if 7* < 7 < 1. 

2. The mathematical model. It is assumed that the dermal tissue is in a circular domain 
{(r, 9; r < L)} and the open wound at time t is a disc {(r, 9; r < R(t))} with initial radius 
R(0) < L. The partially healed tissue is the annulus fi(t) = {(r, 6; R(t) < r < L)}. We 
introduce the following variables: 

• Chemicals: 

w(r, t): concentration of tissue oxygen 

e(r, t): concentration of Vascular Endothelial Growth Factor (VEGF) 
p(r, t): concentration of Platelet Derived Growth Factor (PDGF) 

• Cells, blood vessels and matrix 

m(r, t): density of macrophages 

f(r, t): density of fibroblasts 

n(r, t): density of capillary tips 

b(r, t): density of capillary sprouts 

p(r, i): density of the ECM 

v(r, t): velocity of the ECM 
In homeostasis w = wo, to = mo, / = fo,b = bo and p = po. In the remainder of this paper 
these variables have already been scaled so that wo = too = po = bo = po = 1. 
The continuity equation for the matrix density p is 

^ + V-(pv) = G p (f,w,p), 

where G p (f, w, p) is a growth and decay term of the ECM due to collagen secretion by fi- 
broblasts and degradation by matrix metalloproteinases (MMPs). The specific form of G p 
incorporates the fact that collagen production and maturation require the availability of oxy- 
gen (13 El HU ED , 

k p w p 

W + K wp p m 

where p m is the maximum matrix volume fraction permitted in the partially healed region, 

Prn > 1. 

The partially healed tissue is modeled as a quasi-static upper convected Maxwell fluid 
with velocity v, deviatoric stress tensor given by r = ?y(Vv + Vv T ), where r\ is the shear 
viscosity, and pressure P. The pressure P is generally a function of the matrix density p, and 
is assumed to have the form 

\ 8(p-l), p>l 
0, p< 1. 



The total stress a = r — PI appears only in the boundary conditions. By further assuming 
radially symmetric flow, i.e., v = v(r, t)e r , the continuity equation becomes 

d P + ±JL(r P v) = ^—f(l--t-)-\ pP , R(t)<r<L, (2.2) 
dt rdr K w + K wp p m 



and the non-dimensionalized momentum equation for the ECM becomes (see B28I , supporting 
information) 

1 d fdv\ v _ dP(p) R ( t \ <r<L (23) 
r dr \ dr J r 2 dr ' 

To simplify the analysis and simulations we wish to have a PDE system in which all vari- 
ables are radially symmetric. In order to implement ischemic conditions in radially symmetric 
form we assume that small arcs of length S are cut off from the healthy tissue at r = L and 
that the distance between two adjacent 5 arcs is e. If 5, e — > in such a way that e ~ e" c / s 
where c is a positive constant, then, for any diffusion process with boundary conditions 

— = on the o-arcs, 

dr 

u = g on the remaining arcs, 
the limiting "homogenized" boundary condition is 

du 

(1 - 7)(u - g) + 7— = on r = L 
dr 

for some constant 7 E [0, 1] which depends only on c; 7 = corresponds to healthy tissue 
(i.e., no excision of (5-arcs) and 7 near 1 corresponds to extreme ischemia. 
The equations for the concentrations of oxygen, PDGF and VEGF are: 

dw Id, \ 1 d ( dw\ ,„ JN 

■77- — I ;r— (rwii) = — — I vD w — — ) (2.4) 



dt r dr r dr \ dr J 

+k w b((l - 7)w h - w) 



(X w ff + \wm,m) [1 + ^ WWP 

' \ 1+p 



+ K 



w, 



dp Id, \ 1 d ( dp\ , . . X v ffp 

dt + rdr^ = rdr [ rD >£ ) + " ~l+p ' (2 ' 5) 

<9e 1 d 19/ <9e\ 

— + --x-( re ' 1 ') = rD eTT + k e mG e (w) ~ (A e „n + A eb 6 + A e )e, (2.6) 

dt r dr ' r dr \ dr J 

The equations for macrophages, fibroblasts, capillary tips and capillary sprouts include 
diffusion, generation and death of cells, and chemotactic migration of cells: 

^l + l^( rmv ) = ( rD dm ) 1 9 I r XmpmH(l - m/m m )dp/dr 1 

^ rdr K rdr\ m dr J r dr ^ ^ + ~ " 

+ ^-X m m(l + X d D(w)), 
l+p 

% + 1 * {r/ „) . I » U «) _ 1 » LM-f/uv»/» 1 (2 . S) 

at r dr ' r dr \ dr J r dr \ /, , i a ..,o ,2 1 



\Jl + k sg \dp/dr 
+k f G f (w)f (l - £) - A//(l + ArfD(u>)), 



+ (k nb b + k n n) 



1 



(X nb b + X nn n)n, 



db Id 
dt r dr 



(rbv) 



ld_ f db 

r dr \ b dr 



1 d 



ld_ 

r dr 



AD n bdn/dr 



(2.10) 



r dr 



'l + k sg \dn/dr\ 

AxnbpnH(l - n/n m )de/d r 
t \de/dr\ 2 



-k b G b (w)b(l - b) + G b {w)(X nb b + X nn n)n. 



where the two terms with A (in (12.1 Oi l) represent the fact that sprouts follow tips, and the 
oxygen-dependent functions G's and D are given by 



' 3w, 
2-w, 



< to < 0.5 
0.5 < w < 1 



G p (w) = < 



1 2 

+ l<w<4 

2, to > 4 



'2w, 0<w < 0.5, 

2-2to, 0.5<w<l, 

G e (w) = \ 1 1 

-w - -, 1 < w < 4, 



to > 4 



+ l)w 



K. 



G b = 



(K wp + l)w 



wf 



K, 



wp 



W 



D(w) = 1 - H(5w - 1)H(1 - to/3). 



Here H is an approximated Heaviside function 




Note that in Equation (|2.4[) the supply of oxygen from the vasculature is reduced to 
k w b((l — 7)u>b — w) due to the ischemic condition. The functions G p (w) and G e (w) are 
constructed to reflect the biological effect of oxygenation: moderate hypoxia and hyperoxia 
increase the production of PDGF and VEGF compared to normoxia. Equations (|2.7[) - (|2.9[) 
include chemotaxis flux terms that describe the chemotactic movement of macrophages, fi- 
broblasts and capillary tips. The two terms with A in Equation (|2.10p represent the fact that 
capillary sprouts are dragged along capillary tips. Although the forms of the G functions 
and D function are suggested by biological experiments, our mathematical analysis will not 
depend on the special form of these functions. 

The free boundary r = R(t) is moving with velocity v: 



R(t) = v(R(t),t) 



The boundary conditions at r = L are 

v = 0, 



dw 



(1- 7 )( W -1) +7 L— = 0, 



(2.11) 

(2.12) 
(2.13) 



3v 3c 
(l- 7 ) p + 1 L-£ = 0, {\- 1 )e + 1 L— = 0, 



(1 — 7)m + 7-L 



c*to % m pmH(\ — m/m, rn )dp/dr 
~dr~~~D r , 



^l + k sg \dp/dr\ 



(1 _ 7)(/ _ 1) + 7L| |f xfPfH(i-f/f m) ap/ar 



dr D 



f y/l + k ag \dp/dr\' 



0. 



, <9n Xn pnH(l - n/n m )de/dr 



<9r D„ 



96 



= 0, 



(l- 7 )(6-l)+ 7 Z|-+ 

Jl + k sg \dn/dr 



1 + fc sg |0e/0rf 

AD n bdn/dr Ax n bpnH{\ — n/n m )de/dr 



yl + fcsg |<9e/<9r|^ 



(2.14) 
(2.15) 

(2.16) 

(2.17) 

-0, 
(2.18) 



and the boundary conditions at r = are 



— = P 

<9w; <9e dn db 
dr dr dr dr 



_ dp _ k pb R(t) 
dr D p R 

dm pmH(l — m/m rn )dp/dr 

Or L , , la , Q |2 



^l + k sg \dp/dr\ 



df pfHQ. - f/f m )d P /dr 

f~dr + Xf / — = °' 

y 1 + k sg \dp/dr\ 



(2.19) 
(2.20) 
(2.21) 

(2.22) 
(2.23) 



Equation (|2.21[) represents the fact that secretion of platelets decreases with healing (i.e., as 
R(t) decreases). The initial conditions for R < r < L take the form 



i?(0) = i? , v = 0, p = f = 1, w = l, b = g 
e = m = n = 0, p = Po(r), 



r-Rp 



(2.24) 



where 



f0, 



i 4 1 

9{z) = < ~z 



6' 

u 



z < 0, 



0<*<-, 

1 3 

4 " 4' 
3 

4<*<1, 
z > 1. 



and po(r) has three continuous derivatives and satisfies the boundary conditions ( 12.14t and 
d%2fl . and 

(p' (r)<0 if R <r<R + e , 
\po(r)=0 if Ra + e <r<L 

where < So < L — R . 

In a healthy tissue there is no net growth of ECM, i.e., G p (f, w, p) = 0iff = w = p = 
1, which means that 

x p = ttit- i 1 -— )■ (2 - 26) 

1 + K wp \ p m J 



Similarly 

= X - 

Wb - 1 
A/ 



k w = A " / + A " m , (2.27) 



3. R(t) is monotonically decreasing. Set 

Q(t)= [ yP(y,t)dy, (3.1) 

jR(t) 

where P{r,t) = P(p(r,t)). 

THEOREM 3.1. For any solution ofM.2\- $2?2M there holds: 

R(t) < 0; R(t) < if and only if Q{t) > 0; (3.2) 

R(0)e-T* f <> Q{T)dT < R{t) < i?(0)e~^ f « Q{T)dT . (3.3) 

Proof. Equation ( 12.3b can be rewritten as 



V r V /v\ 
\r J r 



2 



r r 

Integrating over [R(t), r], we obtain 

v r (r,t) - VrW)) + ^ " V ~^- = P(r,t) - P(R(t),t). 
From d2.1 U and ( 12.191 > we obtain 

« r (r,*) + -^i-^ = P(r,t), 



hence 



(rv) r -r^ = rP(r,t). (3.4) 

6 



rv(r,t) — ~^r= I yP(y,t)dy. (3.5) 



Integrating this equation over [r, L] and using ( |2.12| i, we obtain 

- r 2 R(t) 
~2 R(t) 

In particular, at r = R(t), 

i ' , L 2 ~ R(t) 2 R(t) f L , N 
-R(t)R(t) -LL -Al = / y P (y } t )dy, 

or 

W) = -vrm? m - t3 - 6) 

The assertion ( 13. 2b now follows immediately from (13.6b . From (13.6b we also obtain 

from which we deduce the estimate (13.3b . □ 

If we substitute R/R from ( 13.6b into ( 13.41 ) we obtain, after dividing by r, 



this equation will be needed in the remainder of this paper. If we substitute R/R from ( 13.6 
into ( 13.51 ). and divide by r, we obtain an expression for v, 



w(r ' <} = " ^ UTRW Q{t) ~ ' yP{Vl t)dy 



or 



if L 2 -r 2 f r r 2 + Rfi) 2 /" L I 

" <r - '> " ; { UTW J vPly - t)dy - / " pfc "*} (3 ' 9) 

COROLLARY 3.2. Equation ( 12.3b for u together with the boundary conditions A2.12\l , 
\2.19\ and the initial condition v = can be equivalently replaced by the formula ( 13.9b . 
In the remainder of this paper we shall often work with the representation ( 13.9b for v. 

4. A priori estimates. In this section we assume that there exists a classical solution to 
( 12.21 ) - (12.28b for < t < T, and derive a priori estimates which depend on T, but remain 
uniformly bounded for any finite T. We set 

fl T = 9, t) | R(t) < r < L, < 6 < 2tt, < t < T} 

and introduce the following notation: 

Cl~\ a ' 1+a ^ 2 (VLt) is the space of functions u(r, t) with u, D 2 u, D t u uniformly Holder 
continuous in &t, with exponents a in r and a/2 in t; the norm in this space is defined by 

|M| c 2+ q , 1W 2 (At) = IM| l ~(q t ) + \\ D rU\\c^\n T ) + ll- D *«|| c -.-/2 (fiT) 
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where 



.. „ [u(r, i) — u(r',t')l 

Hlc?i« / 2 ( n x ) = IMU-PM + / w su P^ _ \ r _ r ,\ a + lt _ t/ \ a /2 - 

r,t v ' r,t), (r'tf )SO T l r r I T \ l l \' 



(r,*),(r',f ))en 

Similarly we define the spaces (f2 T ), C 1+Q [0, T], etc. 

In the remainder of this paper we shall use the following comparison principle ||7l [T5l . 
Lemma 4. 1 . Let v%, v 2 satisfy 

-^-DAv 1 +g(x,t,v 1 ,Vv 1 )>-^--DAv 2 +g(x,t,V2,'Vv2) in Vt T . (4.1) 

//' 

Q 

Hig^(vi ~ v 2 ) + ll 2 ( Vi - v 2 ) >0 on dn T n{o<t<T}, 

(vi - v 2 )\ t= a > in fi 

where v is the outward normal and fjti, p 2 are nonnegative functions satisfying, at each point, 
either \x\ > or p,\ = 0,/i 2 > 0, then v\ > v 2 in fix- Furthermore, if strict inequalities hold 
in both (|4.ip and (|4.2[) , f/ien «x > ^2 in S^t- 



Lemma 4.2. For any solution of (|Z2| - (|2.2g|> . 

the components w , e, p, m, /, n, 6, and p are nonnegative functions. (4.3) 

Proof. For any small <5 > 0, let us add <5 on the right-hand side of each of the equations 
(|2!4|) - (|2l0|) and each of the boundary conditions (j27T3|) - pTT8]) . (|23I]) - (|2!23]) . replace by 
—d in pT20|) . and increase the initial data of 6, e, m, n,p by 5. We refer to this new system 
as the "J-problem" and to its solution as the "<5-solution". By continuity, each component of 
the (5-solution is strictly positive in Vl to for some to > 0. We claim that all the components 
are strictly positive in ft? for all T > 0. Indeed, otherwise there is a smallest T such that at 
least one component of the (5-solution, denoted by z, vanishes at some point (f, T). We can 
then apply the second part of Lemma |4~T1 with Vi = z, v 2 = to conclude that z(f , T) > 0, 
which is a contradiction. 

The local existence and uniqueness proof given in Sections 4-6 is valid also for the 5- 
problem. The estimates derived there are uniform in S so that, as 5 — > 0, the (5-solution 
converges to the original solution. Hence each component of the original solution is non- 
negative in a small time interval, say < t < t*. We can now repeat the process for t > £*, 
and conclude, step-by-step that each component of the solution is non-negative in fix for any 
T > 0. □ 

LEMMA 4.3. If initially p(r, 0) < p m forR(0) < r < L, then, 

p< p m in fi T . (4.4) 



Proof. If the assertion (14.4b is not true, then there exists a t* > such that p(r, t) < p m in 
Q,f, and p(r*,t*) — p m for some R(t*) < r* < L. Then, along the characteristic curve 
with velocity v, through (r* ,t*), 



(4.5) 



where D/Dt = d/dt + v{d/dr). On the other hand, from J2.2| i and d3.81 l we get, 



Dp 
Dt 



(r«,t«) 

2 „2 



-Ki^.r)-[p{r\t*)-j^-^Q{e)),,u-'.r). 



Since Q(t*) < L ~ R max r P(r, i*) = L P(r*,t*), we obtain 



Dp 
Dt 



(r*,t*) + -fl 



which is a contradiction to (14.5b . □ 

Recall that we have assumed p m > 1. 
LEMMA 4.4. TTzere /ioMs: 



b(V til 

LA±J1 < P(p m ~ 1), \v r (r,t)\ < 2f3(p m - 1), in Q T . (4.6) 



Proof. From Lemma |431 we obtain 

r 2 - P(Tl 2 

yP(y,t)dy<p(p m -l) ^L, 

1 

L L 2_ r 2 



yP(y,t)dy < p( Pm - 1) 
Using these estimates in ( 13. 9t we get 

Substituting this inequality into ( 13.8b and estimating P and Q by Lemma l4~3l we also obtain 

K(r,i)| <2/3(p m -l). □ 



Lemma 4.5. Setting 



knb K + f3[p m - 1] 



N = max <^ - — , , n m ( 

there holds: 

< n(r,t) < N in T . (4.7) 

Proof. We write Equation J2.9t for n in the form 

.Sf[n] = JSf [n] + ^[n] = 0, 

where 

„ . 90 1 9 / 5<A 1 a / x»P0g(l - ^/n m )de/ dr 

^o[<p] = -57 - ~ +«^r + --5- r , = 

at r Or \ or I r or \ , , , , n /a ,2 



^l + k sg \de/drf 



and 

&[4>] = b (\ nb <t> - k nb -—\ + (\ nn <t> + - k„-^—) 

\ 1 + ej \ r 1 + eJ 

By (EH) and Lemmal431 



\(rv) r > -I3[p m -1], 



so that, by definition of N, 



and 



[TV ) 6 

XnnN + - k n —— > X nn N - fi[p m - 1] - k n > 0, 

r 1 + e 



X nb N - k nb — — > X nb N - k nb > 0. 



Since, by ( 14.3b . 6 > 0, we conclude that ,^[N] > and hence N is a supersolution, i.e., 
■Sf(N) > 0. Using also the boundary conditions (12. 17b and (12.20b we deduce, by the com- 
parison lemma, that n(r, t) < N. □ 

LEMMA 4.6. For any T > 0, there exists a constant Ct such that 

< b{r,t) < C T in Q T - (4.8) 



Proof. By the comparison principle, 

< b(r,t) < 6i(r,t) 

where bi (r, t) is a solution of the same equation as b(r,t) but without the quadratic term 
—k b G b {w)b 2 and with the same boundary and initial conditions as for b(r, t). We can write 
the equation for b\ in the form 



dh d 
dt dr 



(rDb^-J +oi (r, t)6i (r, t)+a 2 (r, t) ^ (r, t) + -^- (a 3 (r, t)6i (r, *)) = <u(r, t) , 

(4.9) 

where, by using (14.7b . we find that oi, 02, 03, 04 are all uniformly bounded. From the Nash- 
Moser estimate [ 15] we deduce that, for any < t\ < T, 

\\h\\c^ /HQtl ) <Ct + Or||6i|| £ - {0ti)j (4.10) 

and by interpolation, 

INU~(n tl) < ||6i(-,0)||x,~ +t? /2 (l+ sup \R(t)\ a/2 )\\h\\c^ 



0<r<ti 



< 



||6i(-,0)|| £ - + C*t1 /2 (CT + CtIMl^^)) 



Choosing t\ such that 

C *C T t? /2 = i 
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we obtain the estimate 

||&i|U°°(f3 tl ) < c - 

Repeating this procedure step-by-step, the assertion ( 14. 8t follows. □ 

The above proof can be applied successively to m, f, p, e and w to establish the following 
estimates. 

LEMMA 4.7. For any T > 0, there exists a positive constant Ct such that in VLt, 

< m(r, t) <C T , < f(r, t) <C T , < p(r, t) < C T , 

<e(r,t) <C T , <w{r,t) <C T - (4.11) 

Since b is bounded (by Ct) in Q,t, we can write the equation (12. lot for b in the same 
form as Equation J4.9t for &i and thus derive, by the Nash-Moser estimate, a Holder bound 

l|fr||c=.°/ 2 (o T ) ^ Ct- 

The same bound can be derived for the components n, m, f, p, e and then also for w. Hence, 
we obtain 

LEMMA 4.8. For any T > there exists a positive constant Ct such that 

\\w,p,e,m,f,n,b\\ ca ,a./2^ T) <C T - (4.12) 



Rewriting Equation (12.2b in the form 

Pt + v P r = -^/(l - — ) - A P P - { -^P = Hr, s), (4.13) 

w + K wp p m r 

we proceed to establish a Holder estimate for the function p. 

LEMMA 4.9. For any T > there exists a constant Ct such that 

\\p\\c^{h T ) < Ct- (4.14) 

Proof. We introduce the characteristic curves X, for ( 14.13b . by 

dX r (r, t, s) 



ds 

X r (r,t,t) = 1 



v r (X(r,t,s),s)X r (r,t,s), Vs G [Q,i] 



Using Lemma l4~4l we find that 

|^r(r,i, S )| < e 2/3(p m -D(*^). 

Let J(r, i, s) = p(X(r, t, s), s), so that 

dJ(r,t,s) 



ds 

J(r,t,t) = p(r,t) 
11 



^(X(r,t,«),«), 



Then 



\p(n,t) - p{r 2 ,t)\ 



\ri - r 2 \ 



< 



1 



\n - r 2 \ c 



F{X{r\,t, s), s) - T(X(r 2 , t, s),s)ds 



+ 



\p(X(r u t,0),0)-p(X(r 2 ,t,0),0)\ 
\ri-r 2 \ a 



By the initial condition p(r, 0) = 1 the last term vanishes, and 



F{X{n,t, s), s) - T(X(r 2 ,t, s), s)ds 



|n -r 2 \ a 



< 



< 



< 



J-(X(n, t, s), s) - T(X(r 2 ,t, s), s) ( \X( n ,t, s) - X(r 2 ,t, s)\\ a 



I \X(n,t,s)-X(r 2 ,t,s))\" V ki-ral 
(p 2p(p m -i)(t-s))a f \F(X(r x ,t, s),s)- T{X{r 2l t, s), s)\ 
[ ] Jo \X(r u t,s)-X(r 2 ,t,s))\<* 

Ct f [p(;s)] c? + H;s)] c? + [f(;s)] c? ds. 
Jo 



ds 



Hence 



\ P (n,t) - P (r 2 ,t)\ <Ct + Ct '■' 



Taking supremum over n, r 2 G [R(t), L\, r x ^ r 2 , we obtain 



/ \p(;8)]c?ds. 
Jo 



[p(;t)] C? <C T + C T [ [p(;s)] C ?ds, 

Jo 



and by Gronwall's inequality, 

[p(.,t)]c ? <C T . 
Next, taking t 2 > t\ > 0, we can write 



(4.15) 



/•*2 

p(r,t 2 ) - p(r,*i) = / J 7 (X(r,t 2 ,s),s)ds + p(X(r,t 2 ,h),t 1 ) - p(r,h), 
Ju 



so that 



Since 



p(r, t 2 ) - p(r, h) < C\t 2 - h | + \p(-, h)] c ? \X(r, t 2 , h) - r | c 



\X(r,t 2 ,h) -r\ = \X(r,t 2 ,h) - X{r,t 2 ,t 2 )\ < 



dX 



ds 



1*2 -*l|, 



we obtain 



\p{r,t 2 ) - p{rM)\ <C T \t 2 -h\ a . 
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Combining this inequality with ( 14.151 ), the assertion J4.14I ) follows. □ 

LEMMA 4.10. For any T > there exists a constant Ct such that 

M\c:f(n T ) + \\v r \\c?f(n T ) < Ct- (4.16) 

Proof. The proof follows from the representations of v(r, t) and v r (r, t) in ( 13. 9b and ( 13.8b by 
using Lemma l4~9l and the boundedness of R (from ( 13.3b ). □ 

LEMMA 4.11. For any T > there exists a constant Ct such that 

\\R\\ci+«([o,t]) < C T - (4.17) 

Proof. This follows from (12. lit and Lemma 14.101 □ 

LEMMA 4.12. For any T > there exists a constant Ct such that 
(i) 



\p\ 




,1+c 


• /2 (n r ) 


<C T , 


\e\\ 






/2 (n T ) 


<C t , 


w 




>,1+ 




t <c T 



(i>) 



m\\ c 2+ 


a,l+ 


a/2 (^T 


} <c T , 


f\\ clT 


;,l + o 


/2 (^t) 


<C T , 


n\\ cl + c 


,,l+o 


/2 (O r ) 


<C t , 


b c 2+o 


,l + o 


/2 (A T ) 


<C t ; 



( in) 

\\P\\c 2 r + a - 1+a/2 (Q T ) ~ Ct > 

\\v\\ c 2+ a ,i +a/ 2^ T) < C T - 

Proof. Indeed, (i) follows from Lemmas 14.81 - RTTTl and the parabolic Schauder estimates 
0Q3]. The assertion (ii) follows by the Schauder estimates and (i). To prove (iii) we first 
formally differentiate ( 14.13b in r and apply the proof of Lemma |4~9l making use of Lemma 
14. 101 and (ii). We thus obtain the bound 

lkll CrT /s ( a T) < Ct- (4.18) 

In order to rigorously prove (|4.18|) . we consider the solution p r of the differentiated 
equation (14.13b and derive the estimate ( 14.18b . By integration of the equation of p r with 
respect to r, one can verify that f r p r dr coincides with p; hence dp/dr = p r and ( 14.18b 
follows. 
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Differentiating d3 ,8b in r and using ( 14.181 ) we deduce that 



W V rr\\ C ^/ 2 (n T ) ^ Ct - 

and this allows us to differentiate the equation for p r once more in r. Proceeding as before it 
is then easy to complete the proof of (iii). 

5. Transformation to a fixed domain. In order to prove existence and uniqueness of a 
solution of (12. 2t - (12.28b for a small time interval < t < T, it is convenient to transform 
the system with the free boundary r = R(t) into a system with a fixed boundary, using the 
mapping 

t=T^Wy {r=(l-QR(t)+ZL). (5.1) 

In the new system £ varies in the interval < £ < 1, and for any function u(r, t) = t), 
du 1 du 



dr L- R(t) <9£ ' 



(5.2) 



and 



d ( du\ 1 d ( , du \ . 

d~r Vfr) = (L-R(t))> d! [ r ®di) ' (53) 



du du dudt; du R(t) , 
~dt = ~dt + ~dH~dt = ~dt + L- R(t) ^ ~ 'di' 



(l-0(i-i2(*)) 



1 u. 



Using these formulas we compute 



du Id, \ du 

+ -—s ruv ) = - d -i+ B ^ 



dt r dr 



where 



or, 



where 



B = 



R(t) (t^M 

\> ) at 



1 d 



L - R(t) ^ ' di [L — R(t))r d£ 



(ruv) , 



R(t) 1 d ( , \ 1 d , _ s „. 



L - R(t) 



ld_ 

r <9£ 



(ru(R(t)(Z 



D + v) 



Ku. 



g=gK)= -j|_r ('-o^-«w) _ 
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(5.4) 



Hence 



du Id 



at 



dr 



[ruvi 



du 
~dt 



L - R(t) 



ld_ 

r 9£ 



Ku. (5.5) 



Using (15.2b . (15.3b and ( 15.5b . we can transform the PDEs in Section|2]into the following system 
of equations, where we have, for simplicity, dropped the tilda "~" from all the variables: 



dp 1 d 

m + ^)dl 
i 



k n w , p , 
to + K wp p m 



1 9 / dv 



1 



<9P 



(£ - R{t)Y r(0 d£ 



+ fc w 6((l - 7)w b - w) 
dp 1 9 / , „ % . ,\ 19 



dU r 2 (C) L-R{t)dt 

« ■ rs|( r «H = 4!( rm,(t) ! 

+ k p mG p (w) — ^ - A p p - Kp, 
+ k e mG e {w) — (A e „n + A e j>6 + A e )e — ife, 



+ A, 



dc 



1 9 



<9i r(£)9£ 



(5.6) 
(5.7) 
(5.8) 

(5.9) 

(5.10) 
(5.11) 



15 



dm 1 9 



dt r(£)9£ 



(rg)mM) (5.12) 



1 9 / n dm\ 19/ Xm{t)p m H(l — m/m m )dp/d^ 



( ,,n „ . . om\ 1 a / 



m f7i (1 + \dD(w)) — Km, 



df 1 9 



m ■ rmeM ),M ) (5 " 13) 



\/i + k sg (t)\d P /dt\ 2 

+ k f G f (w)f (i - - A //(l + AdDH) - Kf, 
°" ' (5.14) 



dt r(Z)d£ 

1 9 / 9n\ 19/ Xf(t)pnH(l - n/n m )de/d£ 



r(0 9£ 



yjl + k sg (t)\de/d£\ 2 



+ (A;„&6 + k n n) — (A„ 6 6 + X nn n)n - Kn, 

1 + e 



96 1 9 



9i r(£)9£ 



r(06Af) (5.15) 



1 9 / 96\ 1 9 / AD n {t)bdn/d£, 



r(0 9£ V 9^ r(0 9£ ^ ^7^^ 

1 9 ( AxnWpnHQ-n/n^de/dZ 



r{0 ^ V Jl + k ag (t)\de/d£\ 2 

+ k b G b (w)b(l - b) + G b (w)(\ nb b + \ nn n)n - Kb, 



where 



« 9 



2 ■ 



L-iJ(t) ' S9W (L — R{t)) 

D u {t) = J " for u = w,p,e,m,f,n,b, 
{L - R(t)) z 

Xu{t) = ( L _ fl( t ))2 ' hi u = m,f,n,b. 



The free boundary condition remains as before, namely, 

R(t) = v(R(t),t). (5.16) 
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The boundary conditions at the fixed boundary £ = 1 are 

u = 0, (5.17) 

t 1 -'«"-"+ At =°' <"*> 

(1 -* + r3Swl = °- (5 - 19) 

(1 - 7)(/ - 1) + -2* ( gLPW-///^/^ = o, (5.22) 

(1 - 7)n + f g - ^^-"/""H = 0l (5.23) 

(l- 7 )(6-l) (5.24) 

7 L / 96 AD n bdn/d£ A Xn bpnH{l - n/n m )de/d£ 

L -m\dt ^ 1 + ksg{t) \ dn/d tf ^l + k sg {t)\de/dtf 

and at the free boundary £ = they are 

dv 



dp k p i,R 



R(0) = R ,v = 0,p = f = l,w = l,b = g 
e = m = n = 0, p(£, 0) = p ((1 - £)Ro + £&) • 



0. 



. ;< (L-R(t))P, (5.25) 

<9ui 9e 9n 96 „ 

0, (5.26) 



<9£ d£ <9£ <9<e 

(L-R(t)), (5.27) 



9£ £ P i? 

9m pmH(l-m/m m )dp/d£ 

^ Vi + M*)|0p/0fl 2 

- % g +X/ ^W^S = 0. (5.29) 
d « yJl + k. B {t)\dpldit 

The initial conditions take the form 

'£( L-Rq) \ 

so J ' (5.30) 



6. Existence and Uniqueness. In this section we prove the following theorem. 
THEOREM 6.1. There exists a unique solution of \2.2\ - ( 12.281 ) for < t < oo such that, 
for each T > 0, the estimates of Lemma \4.12\ hold. 
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Proof. We first prove existence and uniqueness for a small time interval < t < r. For 
this proof it will be convenient to transform the system (|2.21 i - ( 12.241 ) into the system ( 15.61 ) - 
(15.30b with a fixed boundary. Set 

G = {0<£<1}, G T = {(£,*);£ e G,0 < i < T} for any T > 0, 

and introduce the Banach space 

Y = {(R(t),p(£, t)); R(0) = Ro, p(£, 0) = 1 with norm 
\\(R,p)\\ Y = \\R\\c 1 +*/*[o,t] + \\(p,Pi)\\c^^(G T )} 

and the ball 

Y B = {(R,p)eY;\\(R,p)\\ Y <B} 

for any B > 1 + Rq. 

For any (R, p) G Yg we wish to solve the system ( 15.71 ) - ( 15.15b with the correspond- 
ing boundary and initial conditions from (15.17b - ( 15.301 ). Denoting this solution by u = 

(w,p, e, m, /, n, b, v) we shall then define (R, p) by 

^-R(t)=v(R(t),t), R(0) = R o , (6.1) 
at 

f + ^-^( r tt)pM)=^r^f(l- A )-\ P p-Kp, p&0) = l, (6.2) 
r(g c*4 V 7 w + K wp Pm 

where 

£l = (dR/dt)(Z-l) + v k = dR/dt ( (1-Q(L-R(t j) _ \ 
L - R(t) ' L - R(t) \ r J ' 

and set 

{R,p) = W(R,p) 

We aim to prove that the mapping W is a contraction mapping, and thus has a unique fixed 
point. 

As in [9 1 one can prove, by a fixed point argument, that there exists a unique solution u 
for < t < t, for t small, and that 

||i4|| c 2+ a ,i+a/2^^ < C, u — (w,p, e, m, /, n, b, v). (6.3) 

The estimate ( 16.31 ) can also be established by the argument used in the proof of Lemma |4.12| 
From d6.1| ) and ( 16.3b we get 

\\j t R\\c^[o,r] < C, (6.4) 

so that 

\\{M,k)\\ c ^ +a/2{0T) <c. 

18 



We next consider ( 16.21 ), and use the same arguments as in the proofs of Lemma |4~9l and |4~T2l 
(iii), to derive the estimate 

\\p\\ c 2+ a ,i+ a /2^^ < C. (6.5) 

From ( 16.4b , ( 16.5b we deduce that 

J \\R\\&+«[o,t]<Ro + Ct, 

{||(p>l)|| crf / 2(eT) <l + ^. (6 ' 6) 

Hence if r is sufficiently small then W maps Yg into itself. 

We next prove that W is a contraction in Yb- Let (i?i, pi) and i?2, P2 be any elements in 
Yb and denote the corresponding solution by Ui = (wi,pi, ei,rrii, fi,rii,bi,Vi) for i = 1,2. 
Set 

(Ri,Pi) = W(Ri,pi). 

As in [9 1 one can show that 

IK - u 2 \\ c 2+ a ,i+ a /2 {G ^ < C\\(Ri - Ri,pi - P2)\\y, (6.7) 
from which one can easily deduce that 

||^(B 1 -^ 3 )|| oa+ - [D)T] <C||(iJ 1 -iia,pi- / oa)||y, (6.8) 

and 

||(M X - M 2 ,Ki - K 2 )\\ c 2+ a , 1+a/ 2 {OT) < C\\(Ri - R 2 ,Pi - P2)\\y- 

Using arguments as in the proof of Lemma |49l and l4~T2l (iii) and noting that p\ — p% = at 
t = 0, we derive the estimate 

\\p~i - p2\\ G 2+a,i+o,/2^ < C\\(Ri - R2.p1 - Pi)\\y- 

Recalling also d6.81 l and the fact that Ri — R 2 — at t — 0, we deduce, analogously to (16. 6t , 
that 

\\(Ri - R 2 ,pi - p 2 )\\y < Ct 1 ' 2 \\{Ri - R 2 , P i - P 2)\\y. 

Hence if r is sufficiently small then W is a contraction. We have thus established existence 
and uniqueness for a small time interval < t < r. 

In order to prove existence and uniqueness for all t > we suppose that such a global 
solution does not exist and derive a contradiction. Suppose that a unique solution exists for 
< t < T but not for a larger time interval. We then use the a priori estimates of Lemma 
l4.12l combined with local existence and uniqueness to extend the solution to a larger interval 
< t < T + t, which is a contradiction. □ 
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7. Ischemic wounds do not heal. In this section we prove that if the parameter 7 in 
the oxygen equation J2.4t and the boundary conditions ( |2.13t - d2.18t is near 1 then R(t) = 
const. > for all t sufficiently large, that is, ischemic wounds do not heal. 

For any function u(r, t) we introduce the integral 



Iu(t) 



ru(r, t)dr. 



(7.1) 



R{t) 



Using ( 12.111 ), ( 12.121 ) we obtain 



d_ 

dt 



R(t) 



ru(r, t)dr 



R(t) "t 



R(t) 



dt 



L du , 
r—dr 

R(t) "t 



L d 

— {ruv)dr, 



R(t) 



dr 



or 



d T , v f L [du Id, A 

dt IM J R ,/{m + r^ ruv) ) dr - 



(7.2) 



This formula will be used in subsequent lemmas. 

For clarity we shall denote the solution u by u 7 , and consider first the case 7=1. 
LEMMA 7.1. There holds: 



I Wl (t)= / r Wl (r,t)dr <Ce- A ^\ C = I Wl (0). 
JRi(t) 



(7.3) 



Proof. Multiplying Equation ( 12.41 ) by r and integrating over r 6 (i? 7 (i), L), we obtain, 



dt 



dw~, 



dw~, 



Ry{t) 



™ 7 (r, t)dr =LD W -^ (L) - R(t)D w (R(t)) 



+ / n k w bj((l - ~/)w b - to 7 ) 

jR(t) 



(A w // 7 + A u , m m 7 ) ^1 + Y^^j + ^wm Wjjdr 



so that, for 7 = 1, 



dt 



and CL3J) follows. □ 

LEMMA 7.2. There holds: 



rfi(r,t)dr as t — ^ cxd. 
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Proof. Multiplying Equation ( 12.81 ) with 7 = 1 by r and integrating over r 6 (Ri(t), L) we 
obtain, after using the boundary conditions ( 12.161 ) and ( 12.231 ), 

J t I h( t ) = f R t r\k f G f { Wl )h (l- A) -\ f h(l + \ d D( Wl ))}dr 

<CI Wl (t)-X f I fl (t)- 
Recalling (17.3b we deduce 

If! (t) < (Cit + C 2 )e- ™™{\ wm ,\f}t asMoo. □ (7.4) 

LEMMA 7.3. There holds: 

r-L 



Ipi(i) = / rpi(r,t)dr — >• as f — > 00. 



Proof. As in the proof of Lemma r7T2l one can easily derive the inequality 

I PX {t) < (C 1 t + C 2 )e~ min{x ^^" ]t ->0 as i^oo. □ (7.5) 

From the definition of in (|3.1[) and Lemma l731 we obtain: 
LEMMA 7.4. 77zere /loZds: 

Q^t) = I Pl (t) = [ rPi(r,i)dr as t -> 00. 
/&i(t) 

We next prove: 

LEMMA 7.5. There exists a constant C such that 

max Wr,t) < Ce~ Xmmt/2 for all t > 0. 

Rl(t)<r<L 

Proof. For 7 = 1, the oxygen equation can be written in the form 

dwi 19/ dwi\ dwi . . 

-^T 3- rD w —— + v— h = 0, 

at r Or \ or J or 



where 

Si(r,t) 



k w h + (X w ffi + Xw m mi) ^1 + Y^~ J + ' Kwm + ^M 7- '*) 



2Q 1 (t) 



L 2 + R 1 {t) 2 



By Lemma 1X41 there exists a t\, such that, when t > t\, 2Q 1 (t)/(L 2 + R\{t) 2 ) < \ wm /2. 
Hence 

Si(r,t) > \ wm /2, 

and by the comparison lemma, 

wAr,t) < max w x {r, ti)e~ Xmmt/2 . □ 

B,l(t)<r<L 
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LEMMA 7.6. There exists a positive constant F*, F£ > f m , such that 
fi < F* for all R t (t) < r < L, t > 0. 

Proof. From Lemmas [7~4l and [731 it follows that there exists a t\ > such that, for all t > t\, 



2 



rQi(t) + k f G f (wi)(l - < Xf/2. 



Using this in ( 12.8b and setting 

fx = max fi{r,t), 

0<t<ti,Rx(t)<r<L 

we deduce by the comparison lemma that 

fi(r, t) < max{/i,/ m } foralli>ii. □ 



We next improve Lemma 173 
Lemma 7.7. 



Proof. By Lemma 17. 4 



max Pi(r, t) — >0 as t — > oo. 

Rl(i)<r<L 



2 

~ (t) < Ap/2 if t>ti. 



L 2 +i?!(i) 2 
Using also Lemmas l7.5l and l7T6l we obtain 

^-<-^Pi + ^-^ max Wl (r,t)e~ K " mt/2 forallt>*i, 

Dt ~ 2 H K wp fl!(*)<r<L n ' ; ' 

where D/Dt is the derivative along the characteristic curves, and assertion of the lemma 
follows. □ 



Lemma 17771 implies that Pi = for all t sufficiently large, say, for t > T*. Hence also 
Qi{t) = if t > T*. Recalling ( 13.61 ) we conclude: 



Lemma 7.8. There exists R\ > and T* > such that 



Ri(t) > R{ for all 0<t<T*. 
#i(t) = Rl for all t>Tf. 



We next extend this result to all 7 near 1 . 

THEOREM 7.9. For any < 1 - 7 < 1, there exists R* > and T* > such that 



Ry(t) > R* for all < t < T^, 
Ry(t) = R* for all t>T*. 
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Proof. Since the estimates of Lemma 14. 121 hold uniformly in 7, any sequence 7$ — > 1 has a 
subsequence for which the solution u 7 of J2.2t - ( 12.281 ) converges in J7 r , for any r > 0, to a 
solution ui of ( 12.21 ) - ( 12.281 l with 7=1; the convergence is in the norms of Lemma ( 14.12l > 
with a replaced by any < (3 < a. Since (by Theorem 16. U the solution of ( 12.21 i - ( 12.281 ) 
with 7 = 1 is unique, we conclude that as 7 — > 1 the solution u 7 converges to u±. It follows 
that 

fh(r, ti) < u; 7 (r, ii) < 7/0, / 7 (r, ti) < F* + 1, i? 7 (t~i) > R{/2 

if ^ is large enough, provided 7 € (70, 1) and 1 — 70 is small enough; here r/ is chosen small 
enough so that 



2 m k p {Fl + 1) 3 



< tA p , (7.6) 



Let [ii , t 7 ) be the maximal interval such that 

p 7 (r,i)<l, Vte[fi,t 7 ), 

We want to prove that f 7 = +00. Noting that Q 7 (t) = for fi < t < i 7 , we also have 

v 7 (r,t) = and R 7 (t) = i? 7 (ti) forfi < t < t 7 . 

Let W(r,t) = r n (r-R) 2 +-q Q where .R = fl(fi), ^ = (l-7)/Iand^= 2jL(L-R). 
Then (dW/dr) (R, t) = and 

cW 

(l-7)(W-l)+7£— >0 at r = L 
or 

if 1 — 7 is small enough. Also 

W t - £>„AW > - 7 )w 6 - TV) - A™W if m < %, 

that is, if 7 is restricted to a very small subinterval (71, 1) of (70, 1). By the comparison 
lemma we then get 

to 7 (r,f) < W(r,*) for t e [fi,iy), 

and, in particular, 

iu 7 (r,t) < 2^ for < e [fi,t T ). (7.7) 
From (T2~2l . (|7\6|) and (f777|> we then obtain, for 7 € (71, 1), 

J£^ X P(J-P) for t G [fl, try), 

so that 

3 

Pj{r,t) < -, for i e [ii,i 7 ). 

This implies that f 7 = +00, and consequently Q 7 (i) = for all t > t\, and the theorem 
follows. □ 
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8. Wounds that do not heal. A wound may be considered to be (completely) healed if 
R(t) as t — > oo. Indeed, biologically, if R(t) becomes smaller than, say, 10 fim (which 
is roughly the diameter of a cell), no cell can move in to occupy the remaining open space of 
the wound. We say that a wound does not heal if 

lim IL(t) = R* > 0. (8.1) 

In Section [7] we proved that if 7 is near 1 then the wound does not heal and, moreover, 
i? 7 (t) becomes constant for all t large enough. In this section we want to explore some of the 
implications of (|8.ip . In particular we show that in wounds that do not heal, the concentration 
of oxygen and the density of ECM cannot exceed those of a healthy tissue as t — > 00. 
Theorem 8.1. IfA8.1i holds then 

Umsup / 7 (r,t) < f m , (8.2) 
lim sup w 1 (r, t) < max{l, (1 — 7)^}, (8.3) 

t— >oo 

lim ess sup p^(r,t) < 1. (8.4) 

t— >oo 

Proof. By ( 13.61 ) and ( 13.11 ). the function Q(t) satisfies: 

L 2 + R 2 (t) . 

QS) = 2 p (?) ^rW- (8 - 5) 

Integrating over (0, 00) and recalling (18.11 ). we conclude that 

^(0)r2,^2 r.2 / uin\\ R 2 (Vi\ — (R*\ 2 



r^,,, f^'tf + z 2 , L 2 fR(0)\ R 2 (0)-{Rt) 2 



< 00. 



We next prove 

f(r, t) < C, for Rj(t) <r<L, < t < 00. (8.6) 
By d3 ,8b we can rewrite the left-hand side of ( 12.81 ) in the form 

l>J 1 + ft ( P-y(r,t) - r2 2 n ,^ .Q,(t)) ■ (8.7) 



dt 7 dr J1 \ 7 v 5 / L 2 + R 2( t y 

Hence the function g(t) = f m e^° T 1 ^" 1 ^ ds is a supersolution of ( |2.8t and, by the comparison 
lemma, 

/ 7 (r,t) <g(t), R(t)<r<L, t > 0. 



Since, by (18.5b . g(t) is uniformly bounded, ( 18.61 ) follows. 
We next prove that 



? 7 (t)|<C* foralU>0. (8.8) 
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We write (12.2b in the form 



OT ror\ J r or u> 7 + K wp p n 



or 



<9P 9 



By (|3.8|) (or (|4.7[) ). (|8.6p and the bound p 7 < p m , we see that the right-hand side is 
uniformly bounded in (r, t). Hence, by integration over R(t) < r < L, 

r — ' dr is uniformly bounded. (8.9) 

R(t) 9t 

Next, by the definition of Q 1 (t) in (|3.1[) . 

dt 



OP 

Q 7 (t) = -R(t)R{t)P{R{t),t)+ I r l^ dr, 



R(t) 

and hence, upon using (|8.9[) and the uniform boundedness of R(t), the assertion (|8.8p follows. 
From ( 13.6b and ( 18.8b . we obtain the estimate 

|i? 7 (i)| < C. (8.10) 

Using the interpolation estimate (see OH, Page 48) 

and noting that the last factor converges to zero as t* — > 0, we obtain 



lim \\R r \\ o ~[t;r+i]=0, (8.11) 



and then, by d3.61 >. also 



lim ||Q 7 ||o-[t*,f+i] =0, V0<a<l. (8.12) 
From ( 18.12b and (18.5b it easily follows that 

Q 7 (t) -> as t -t oo, (8.13) 
hence there exists a T > such that 

^Q 7 W < Y if * > T - (8-14) 

Writing the left-hand side of (|2.8p in the form (|8 and using (|8.13[) . we can then apply the 
comparison lemma to / 7 to conclude that 

A(r,t)</ m + max / 7 (r, T) e -^*- T ), 

R~,(T)<r<L 

and hence (18.2b follows. 

Similarly one can prove, by comparison, the estimate ( 18.3b . Finally, (18.4b follows from 
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9. Simulations and a conjecture. We simulated the radius R 7 (t) of the wound for 
different values of 7 using the nondimensional parameters of the system (|2.2[) - (|2.28|) that 
were chosen on the basis of experimental results ||28l . In Figure |9TT| we present simulation 
results in the original dimensional variables with L = 7.5 mm and initial wound radius 
Rq = 4 mm. The computation was manually stopped when the wound became 98% closed. 
From the figure we see that as 7 increases, the wound closes slower, and when 7 is close to 
1, the wound radius stops decreasing after a certain time. 

Radius of the wound 

5r 




10 20 30 40 50 60 
Time in days 



FIG. 9.1. The radius of the wound as a function of time for different values of-y. From left to right: 7 = 
0, 0.1, 0.2, . . . , 0.8, 0.9, 0.92, 0.95, 1. Other parameters used are the same as in [281; the nondimensionalized 
values are: L = 5, R = 8/3, p m = 2, K wp = K wf 0.25, k p = 5/16, A p = 0.1, [3 = 10, D w = 0.5, 
D p = D e = 1, D m = D f = 5 x lO- 2 , D n = 10~ 3 , D b = 7 x lO" 4 , X m = Xf = 0.1, x„ = 1, m m = 
f m = n m = 10, A = 0.1, w b = 2, k w = 4.39, \ wf = 0.227, \ wm = 4.16, A d = 2, k f = 5.78 x 10~ 3 , 
Xf = 5.2 x 10- 3 , k nb = k n = k b /10 = 2.16 x 10" 2 , X nn = l00X nb = 2.25, k sg = 6.25 x 10" 2 . 

We conjecture that if the parameters of the system i2.2\ - (12.28b are chosen on the basis 
of experimental results, as in [28], then there exists a parameter value 7, such that ( 18. U holds 

if 7* < 7 < 1 and 

lim RJt) =0 if 7 < 7*. 

If this conjecture is true then, in particular, 

lim Ry(t) =0 if 7 = 0. 



But even this assertion is still an open question. We can only prove, for the system (|2.2[) - 
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(|2.28|) , with general parameters, the following result. 
Theorem 9. 1 . Ifj = 0, then 



p(L,t)>l, 0<t<oo, (9.1) 
R(t) < 0, < t < oo, (9.2) 
Q(i) > 0, < t < oo. (9.3) 



Proof. Using the boundary conditions w(L, t) = 1, f(L, t) — 1, v(L, t) = and (|2.26|) . we 
obtain from (12.21) at r = L the relation 



dt dr 1 + X^.p V PmJ 1 + V P 



1 + K wp 



and from (HI 



Hence, 

= -c (p(L,t) - l) -ML,t)( P (L,t) - l) + + J^g )2 Q(t). (9.4) 

where Co is a positive constant. 
Using the initial conditions 

p(r,0) = l, w(r,0) = l, /(r,0) = l, m(r,0)EE0, u(r,0) = for.Ro<r<L, 
b(r, 0) = 1, p(r,0)=0 for i? + £o < r < L, 

in ()2.2p and (|2.4|l and recalling the relations (12.26b and (12.271 1. we find that 

0p(r,O) 



dt 
dw(r, 0) 



9* 

Using (12.25[) we also obtain (upon recalling (|2.28|) ) that 



0, Rq < r < L, (9.5) 
0, Rq + £ < r < L. (9.6) 



^^ = M m i-j-)p' (r)/^l + k sg (^) 2 >0, R <r<R + e . (9.7) 

Differentiating Equation ( 12.2b in f and using d93l - and the C^ a ' 1+a/2 regularity 
of to, we deduce that 

^^>0 forO<i?o + e -r «l 

This implies that , for < R + e - r < 1 and < t <C 1, 

p(r Q ,t) > 1, 
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and hence 



Q(t)>0 forO<i<l. (9.8) 

Since p(0, L) = 1, from (|9"l4|) and {93J it follows that 

p(0,t)>l (9.9) 

for all < i < oo. This in turn implies that Q(t) > 1 for all < t < oo, hence ii(t) < 
and (by fl9g)) f) > 1 for a// < t < oo. □ 

10. Conclusion. In this paper we established existence and uniqueness of a solution to 
a free boundary problem which models ischemic wound healing. The ischemic condition is 
described in terms of a parameter 7 (0 < 7 < 1) which appears as a coefficient in a Robin 
boundary condition for the various cells and chemical densities. We also proved that under 
extreme ischemic conditions (7 near 1) the open wound stops decreasing in finite time. When 
the parameters of the system are taken on the basis of biological experiments, simulations 
show that there is a parameter 7* such that the wound heals if < 7 < 7* and does not heal 
if 7* < 7 < 1. This assertion remains a challenging mathematical open problem. Future 
work should include the introduction of pressure and diabetic conditions in ischemic wounds, 
as well as inflammatory conditions. 
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